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SYMBOL 
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PSD 


NOMENCLATURE 


DEFINITION (units, U.S. Customary) 


Distance of the rotor disk’s centerline from the inboard 
bearing support. (in) 

Influence coefficient due to an unit load applied at the 
location of the rotor disk. (in/Ibf) 

Distance of the rotor disk's centerline from the outboard 
bearing support. (in) 

Diameter of the rotor shaft. (in) 

Modulus of elasticity of the rotor shaft. (psi) 

Discrete Fourier Transform. (mils sec) 

First moment of inertia of the rotor shaft. (in4) 
Length of the rotor shaft. (in) 

Mass of the rotor disk. (Ibm) 

Mass of the rotor shaft. (lbm) 

Number of scans. 

Power spectral density. (mils/HZ) 

Center of mass offset due to a mass imbalance. (in) 
Whirl orbit major radius. (mils) 

Cartesian coordinates. 


Shape function as a function of axial position. (in) 
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SYMBOL 


W22 


GREEK SYMBOLS 


DEFINITION (units, U.S. Customary) 


Weight density of the rotor shaft. (Ibf/in3) 

Phase lag between the maximum rotor displacement and the 
mass imbalance. (degrees) 

Density of the rotor shaft. (Ibm/in3) 

Rotor rotational speed. (RPM) 

First lateral natural frequency of the rotor assembly. (RPM) 
First lateral natural frequency of the rotor shaft only. 

(RPM) 

First lateral natural frequency of the rotor mass mounted on 


the rotor shaft in the absence of other masses. (RPM) 
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I. INTRODUCTION 


Rotordynamics is a continuously growing field, essential for the design of safe, 
reliable machinery. As the demand for higher power, higher performance machinery 
increases, the need to accurately predict and analyze the dynamic response of rotating shafts 
becomes increasingly important. 

The many uses of rotating shafts includes many industrial, space and military 
applications. Examples of such use are steam and gas propulsion turbines, turbine 
generators, power-transmission shafting, pumps and compressors. These machines 
operate through conditions of varying speeds and loads. As technology advances, these 
machines are becoming lighter in weight, higher in speed, and more power dense 
machines, making the shaft more flexible. Thus, making understanding of the various 
complex rotordynamic phenomena which occur more important. 

Inherent to rotating shafts are numerous rotordynamic phenomena associated with 
radial or lateral, torsional and longitudinal vibrations. The list of rotordynamic phenomena 
is extensive. Several examples include resonances at critical speeds, synchronous shaft 
whirl, chatter and rotordynamic instabilities. The causes of the different rotordynamic 
phenomena easily oumumber the types of phenomena, examples of which are residual rotor 
unbalance, shaft bow, rub conditions, bearing support asymmetry and coupling, and oil 
whirl conditions. It is the accurate analytical modeling of such complex conditions which 
challenge the designer. 

In the past, the most extensive portion of the literature on the dynamics of rotating 
shafts was concerned with determining critical speeds and natural frequencies. It followed 


that most rotating machinery had been designed to operate below the first critical speed 


which usually assured safe, reliable operation. However, most modem rotating machinery 
operate at speeds greater than the first critical speed. Designers now face a larger, more 
complex scope of rotordynamic problems, that is, the prediction of rotor response while 
passing through lower critical speeds, the extent of unstable regions and the associated 
stresses developed. 

The advances in rotordynamic theory, computer simulations and experimental 
equipment are all essential in the accurate prediction of rotordynamic phenomena. As 
analytical methods for accurate prediction of rotordynamic response become increasing 
complex, it becomes apparent that an experimental platform becomes an important 
component for validation of analytical models and direct computer simulation. An 
experimental facility that effectively duplicates rotordynamic phenomena, found in modern 
rotating machinery, is extremely useful for laboratory investigation and educational 
demonstration. 

Therefore, the main objective of this research was to develop an experimental facility 
that effectively duplicates the rotordynamic phenomena associated with rotating machinery. 
The rotordynamic platform consisted of a Nevada Bentley Rotor Kit with a computerized 
data acquisition and analysis process. The use of a self-designed LabVIEW Virtual 
Instrument allows for a flexible data acquisition process. 

Following completion of the experimental facility, it was desired to utilize the facility 
in determining rotordynamic characteristics of the rotor assembly during subcritical, 
resonant and supercritical conditions. Experiments conducted included construction of 
cascade plots for a simply supported, single disk configuration with the rotor disk 
positioned at 1/8, 1//4 and 1/2 axial positions. 

Analytical modeling and prediction of the rotor assembly's first lateral natural 
frequency was conducted while moving the position of the rotor disk along the length of 


the shaft. Rayleigh's energy method and Dunkerley's formula were utilized for the upper 
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and lower theoretical bounds for the frequency respectively. Experimentally measured 
natural frequencies were then compared to this bounded region. 

A Mechanical Engineering vibrations laboratory was designed to investigate the 
phenomena of synchronous whirl due to a mass imbalance. A demonstration of this 
experiment showed the rotor assembly experiencing the unique phenomena of backward 
whirl,as well as, oil whip 

Finally, a summary of the supporting conclusions and recommendations are given. 
This facility can effectively duplicate rotordynamic phenomena and should be utilized for 


further investigation of these phenomena. 


It. EXPERIMENTAL FACILITY GENERAL DESCRIPTION 


The experimental facility consisted of a versatile computerized data acquisition process 
which precisely measured the rotor position and rotor displacement of a Bently Nevada 
Corporation Rotor Kit. The rotor kit is a model of high speed rotating machines that 
effectively duplicates the associated rotordynamic phenomena. Data acquisition was 
performed utilizing an IBM compatible personal computer with an internal data acquisition 
board and an external general purpose termination breadboard. The data acquisition 
process was controlled by a software driven Virtual Instrument (VI) developed using the 
software National Instruments LabVIEW. Proper integration of these highly technical 
components enhanced the development of this experimental facility which serves as a test 
bed for rotordynamic phenomena of rotating machinery. See Figures 1 and 2 for 


photographs of the facility and Figure 3 for a block diagram of the facility. 


A. EXPERIMENTAL FACILITY HARDWARE 
1. Bently Nevada Long Base Rotor Kit (Model 24755) 

The Bentley Nevada Long Base Rotor Kit may be assembled to represent various 
operating configurations and conditions. The basic kit includes a rotor and base assembly 
with motor and rotor disk; external power supply and speed control, and sensing probes 
(eddy current proximeters) used to measure rotor shaft displacement. A brief description of 
the major components used for a simply supported, single disk rotor assembly follows. 


This configuration was used for all experiments conducted within. 
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Figure 1. Photograph of the Experimental Facility. 
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Figure 2. Photograph of the Bently Nevada Rotor Assembly in a Simply 
Supported, Single Disk Configuration. 
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Figure 3. Block Diagram of the Experimental Facility. 


a. Motor Speed Control & Power Supply 

The motor speed control and power supply provides power for both the drive 
motor and the eddy current proximeters. The drive motor's speed can be vaned from 
approximately 200 RPM (slow roll) up to a maximum of 11,600 RPM. The rate and 
direction of speed ramping, either for increasing or decreasing speed. can also be selected. 

b. Rotor and Base Assembly 

The Rotor and Base Assembly (Figure 4) consist of a rigid mounting base, 
both an inboard and an outboard bearing support, drive motor, X-Y probe mount, rotor 
shaft and rotor disk. The bearings are Oilite bushings (oil impregnated. sintered bronze). 
The rotor disk, bearing supports and probe mounts can be positioned anywhere along the 
rotor shaft using predniled holes in the mounting base. 

The rotor shaft is connected inline to the drive shaft motor with a flexible 
coupling. A single notched disk and a twelve notched disk, which provides feedback to the 
motor speed controller, are mounted on the drive motor shaft. These disks allow for rotor 
position monitoring during operation. 

Specific dimensions and material properties of the components of the Rotor 
and Base Assembly are found in Appendix A. 

c. Eddy Current Proximeters 

Rotor shaft displacements are measured by eddy current proximeters. The 
rotor kit has the capability to power four proximeters (not including the !2 notched 
feedback proximeter). Proximeter output, in volts, is a function of separation between the 
proximeter and the rotor shaft. A proximeter calibration curve can be found in Appendix 
B. The proximeters are found to be linear between the range of 2 and 55 mils with a output 


rating of 240 mV/mil. 


anQ PROBES _ 
ae : a | 


20fCR wASS = 


t~-* 2@Q08f wOUNT = YE 





mOLl? ZVENLY 378029 





me os 2 4 TN TIRE | NGik 
ae 2 ee CO 
Ney Ss 
“6 POR 12 ned ea 5° aby 
se = NS a yA o 2 
~ / s X | ~ pe ss. NG 
Re | a“ NACARC 5é42ING RCUS 
\ 


- wOTOR $PeZD 


UOTOR, <7t — ZONTROL 

ZONTROL 7ABLE eu 
LP : ie - 7$ 
nw Note, Rotor <its smown without provective sve! 


Figure 4. Bently Nevada Rotor Assembly Diagram. 
(Reprinted with permission from the Bentley Nevada Corporation) 


The configuration and signal connections for the proximeters are explained in 
Appendix C. This configuration was used for all experiments conducted. 

Though not part of the rotor kit, an electrical device was designed to convert 
selected proximeter output to a TTL signal. The converter can handle two signals. TTL 
signals can be used as triggering and clocking signals for the data acquisition process, as 
discussed later. 

2. Data Acquisition Hardware 
a. NI General-Purpose Termination Breadboard (SC-2070) 

The outputs of the proximeter transducers are connected to a National 
Instrument general-purpose termination breadboard (SC-2070). This 1s a device allows for 
the connection of analog and digital signals, as well as, triggering and scan rate signals to 
the data acquisition board. The SC-2070 is connected to the data acquisition board via a 
50-Pin Input/Output Connector. Silk-screened component locations for resistors and 
capacitors are available to build filters and voltage dividers. The board also has a general 
purpose breadboard area for custom circuit design. [Ref. 1] 

The board has the capability to connect eight differential or 16 single-ended 
analog input channels to the data acquisition board. For all experments, differential 
configuration was utilized, thus limiting input eight to channels. Channel designation is 
found in Table I. 

When required, extemal trigger signals are connected to the board's EXT 
TRIG* terminal and external scan rate or clocking signals are connected to the board's 
OUT2 terminal. The grounding connection for both of these external signals is the digital 


ground or DGND terminal. 


TABLE [. SC-2070 DIFFERENTIAL INPUT 
CHANNEL DESIGNATIONS 


Channel Designations 
Differential Input 
Channel Terminals 
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b. NI Data Acquisition Interface Board (AT-MIO-I6F-5) 

The NI Data Acquisition Interface Board (AT-MIO-16F-5) is an internal high- 
performance mutifunction analog, digital and timing I/O board for the personal computer. 
It uses 12-bit sampling ADC with the capability to sample data at a 200 ksamples/sec over 
the 16 single-ended or eight differential analog input channels described earlier. It has the 
capability for either digital or software triggering. The board is installed in a 16-bit 
expansion slot in the computer. It has an onboard DAC reference voltage of 10 V, called 
AGND, any signal greater than +12 V with respect to this reference voltage will saturate 
internal components and result in distorted data. Configuration and control of the data 
acquisition board is performed utilizing NI LabVIEW for Windows software. See 
Appendix D for actual configuration procedures and default parameters. [Ref. 2] 

c. DataStor IBM Personal Computer and HP Inkjet Printer 

The DataStor IBM personal computer is a 486DX2 computer with CPU 
speeds of 66 MHz or 20 MHz, 16M RAM anda515M hard drive. Connected to the to the 
parallel port LPT1 is a HP Inkjet 550C printer. 
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d. Optical Stroboscope 
The use of an optical stroboscope allowed for the accurate determination of 
rotor speed since the motor speed controller does not have precise speed indications. 


Reflective tape was placed of the face of the rotor disk. 


B. EXPERIMENTAL FACILITY SOFTWARE 
1. Computer Operating Systems 

The disk operating system installed on the computer is MS-DOS 6.0. 

Additionally, Windows 3.1 is installed and required for the data acquisition software. 
2. Data Acquisition Software 

The software used to control the data acquisition process was National 
Instruments LabVIEW for Windows version 3.0. LabVIEW is a program development 
application which uses a graphical programming language, called G, to create programs in 
block diagram form. It is designed specifically for data acquisition and instrument control. 

LabVIEW programs are called virtual instruments (VI's) because their appearance 
and operation imitate actual instruments. There are two main parts of a VI, the front panel 
and the block diagram. 

The front panel is the user interface portion of a VI. It acts like the panel of an 
actual instrument. It contains controls and indicators. Controls simulate instrument input 
devices. They supply data to the block diagram. Indicators simulate instrument output 
devices. They display data acquired or generated by the block diagram. 

The block diagram is the graphical source code or program of a VI. It is 
constructed by ‘wiring’ together operators, functions or even SubVI's, which ultimately 
control the flow of program execution. 

The virtual instrument designed for the control of the data acquisition processes 


for the experiments to follow is called ROTOR.VI. This VI is a hardware timed, software 
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triggered, circular buffered I/O process. This VI was utilized for the acquisition, display 
and storage of proximeter output. Required inputs. operating procedures and VI outputs 


are described in Appendix E. 
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WT. ROTOR CHARACTERISTICS: CASCADE PLOTS 


Cascade plots are effective tools that can be use to determine the general fundamental 
rotor characteristics and the rotordynamic phenomena experienced during subcritical, 
resonant and supercritical operation of the rotor assembly. Cascade plots were constructed 
for the rotor assembly for several positions of the rotor disk. A cascade or waterfall plot is 
constructing by ‘stacking’ the power spectral densities of the time trace signal 
corresponding to rotor speeds over the desired range of operation. Information such as the 
first lateral natural frequencies, excitation of rotor sub-harmonics and super-harmonics, oil 
whirl and oil whip conditions, as well as, many other phenomena can be observed from the 


plots. 


A. ROTOR ASSEMBLY CONFIGURATION 

The rotor kit was configured as a simply supported, single disk configuration, with a 
mass imbalance attached. The inboard and outboard bearing supports were set 17 inches 
apart, leaving one half inch overhang on each end. The rotor disk had a two gram mass 
imbalance attached. Experiments were conduct with the rotor disk was positioned at 8.5, 
4.25 and 2.125 inches from the inboard bearing support. 

The virtual instrument used to control the data acquisition process was ROTOR. VI. 
For an explanation of ROTOR.VI, see Appendix E. The number of scans was 1024 with a 
sampling frequency of 2048 samples/sec. The x and y-displacements versus time data 
were collected for rotor speed ranging from 300 RPM to 10,000 RPM, incrementing every 
200 RPM. Figure 5 is an example of the y-displacement time traces for a rotor speed of 


5560 RPM. 


ee 


The time domain signals of the y-displacements voltages were formatted into a matrix 
with rows corresponding to each rotor speed (run number) and columns corresponding to 
the scan number. The DC component of the signal was subtracted, resulting in zero mean 
data, and then converted from voltage to mils using a calibration factor. Each time domain 
data point was then converted into the frequency domain by taking the Discrete Fourier 
Transform (FFT). Since dealing with FFT's, it is required that the number of scans be a 


multiple of two. The power spectral density (PSD) was then calculated using 


PSD = —— (1) 
where N = number of scans (1024). 

The resulting units of the PSD matrix are mils per Hz. Figure 6 is an example of the 
PSD curve of the y-displacement voltage for a rotor speed of 5560 RPM. Plotting the 
matrix in 3-D by using a MATLAB Waterfall command, plots the PSD curve for each rotor 
speed successively. This results in PSD magnitude along the vertical y-axis, frequency 
along the x-axis and rotor speed along the z-axis. Using a view command of view(0,45), 
rotates the plot such that a 2-D cascade plot is obtained with the frequency along the 
horizontal axis, PSD magnitude along the left vertical axis and rotor speed along the nght 
vertical axis. 

Cascade plots were generated for each rotor mass location, over the operating speed 


range. 


B. DISCUSSION OF TITE RESULTS 
1. Rotor Disk at a=8.5 inches: Centered 
Cascade plots with the rotor disk position at the center of the rotor shaft are shown 


in Figure 7. These plots clearly show a synchronous whirl condition along the first order 
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response. For this configuration, peak response indicates the first lateral natural frequency 
(Wn) to be 2250 RPM. Additionally, the second order (2x) super-harmonic can be seen 
for a majority of the running speeds. For rotor speeds near Wn, the second, third (3x) and 
fourth (4x) order super harmonics are also excited. Also, as the rotor speed approached 
twice Wn, nonlinearities develop causing the one half order (1/2 x) sub-harmonic to be 
excited. This could be the beginning of an oil whirl condition. 
2. Rotor Disk at a=4.25 inches: 1/4 Axial 

The rotor disk's distance from the inboard bearing was decreased in half. The 
cascade plots with the rotor disk positioned at a=4.25 inches are shown in Figure 8. These 
plots show all phenomena as in for the plots for a=8.5 inches described above. For this 
configuration, peak response indicates Wp to be 2700 RPM. The oil whirl condition along 
the (1/2 X) subharmonic continues though indications are not as prevalent. 

3. Rotor Disk at a=2.125 inches: 1/8 Axial 

The distance from the outboard bearing was once again decreased in half. The 
cascade plots with the rotor disk positioned for a =2.125 inches are shown in Figure 9. 
These plots show all phenomena as in for the plots for a=8.5 and a=4.5 inches described 
above. For this configuration, peak response indicates Wy to be 3700 RPM. Additionally, 
an second critical speed was observed along the synchronous speed line at 8400 RPM. 
The amplitude between critical speeds remain relatively high. For speeds greater than the 
second critical speed, the amplitude decreases significantly. Once again, as the rotor speed 
approaches twice Wn, the one half order (1/2 x) subharmonic is excited. At speeds slightly 
greater than Wn, this unstable condition continues along a frequency equal to Wp vice along 


the one half order subharmonic. This follows the characteristics of an oil whip condition. 
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Figure 5. Time Trace for the Y-Displacement for a Rotor Speed of 5560 
RPM. Rotor Assembly in a Simply Supported, Single Disk 
Configuration, Rotor Length of 17.0 Inches, Disk Located 8.5 Inches 
from the Inboard Bearing Support with a 2 Gram Mass Imbalance. 
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Figure 6. The Power Spectral Density Plot of the Y-Displacement Time 
Trace Shown Above in Figure (5). 
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IV. ANALYTICAL PREDICTION OF THE FIRST LATERAL MODE 


The position of the rotor disk and its effect on the first lateral natural frequency was of 
interest. Theoretical predictions for the bounded region of Wn was first completed using 
Rayleigh's energy method for the upper bound and Dunkerley's formula for the lower 
bound. An experimental measured wp was determined for the rotor assembly as the 
position of the rotor disk moved along the rotor shaft. Comparison of the experimental wn 


and the theoretical bounded range was then conducted. 


A. ROTOR ASSEMBLY CONFIGURATION 

The rotor kit was configured as a simply supported, single disk configuration, with a 
mass imbalance attached. The inboard and outboard bearing supports were set 14.5 inches 
apart, with the outboard bearing support flushed with the rotor shaft's end. The rotor disk 
had a two gram mass imbalance attached. 

The rotor disk was positioned at positions symmetric to the center of the rotor shaft 
(a=2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 7.5, 8.5, 9.5, 10/5, 1£-5, 12.5 inches). Thesyramat 
instrument used to control the data acquisition process was ROTOR.VI. The default 


parameters for the VI inputs were utilized. 


B. THEORETICAL UPPER BOUND: RAYLEIGII'S ENERGY METHOD 
Rayleigh's energy method for determining the first lateral natural frequency of a 
system provides a theoretical upper bound for the frequency. Rayleigh’s energy method 


neglects damping and assumes that the total energy of a system is conserved. 


Tmax = Umax (1) 


The rotor assembly was modeled as a simply supported beam (pinned-pinned) with the 


rotor mass treated as an applied load, as seen in Figure 10. 


(x into paper) 





t=a+ bp 


Figure 10. Rotor Assembly Model. 


As derived in [Ref. 3], the maximum kinetic energy and maximum potential energy of 


the shaft and disk are 
[i 
y 2 m 2 
—— 5 | (¥@en) dz +>(y(z =aho,) (2) 
0 
eae 2 
Blo ds 
Uae = sien dz (3) 
0 


where E =the modulus of elasticity of the rotor shaft 


I= the first moment of inertia of the rotor shaft 


a 


L = length of the rotor shaft 
y = weight density of the rotor shaft 
m = mass of the rotor disk mass 
y(z) = shape function 
Substituting equations (2) and (3) into equation (1) and solving for the first lateral 


natural frequency yields 


— 





(4) 
vf y(Z)W _ ) *d2 +m (y(z=a)w,) \e 


i 
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The shape function used was the static deflection equations for a simply supported 


( 
| 
wo. =! 
ee ee 
| 
| % 


beam with a concentrated load 


y(z) = io (i a b” | for z <a (5) 
y(z) = ary (V -b7)z = 3 e(z = P| forz2a (6) 


The first lateral natural frequency was determined while varying the location of the 
rotor disk utilizing a computer program for MathCAD 3.0 for Windows. MathCAD's has 


the ability to evaluate derivatives and integrals directly. 
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C. THEORETICAL LOWER BOUND: DUNKERLEY'S FORMULA 
Dunkerley's formula for determining the first lateral natural frequency provides a 
theoretical lower bound for the frequency. 
Asin the Rayleigh's Energy method, the rotor assembly can be modeled as a simply 
supported beam (pinned-pinned) with the rotor mass treated as an applied load, see Figure 
10. 


As derived in [Ref. 4], Dunkerley's formula for the rotor assembly can be written as 


SG a 1) US ea a (7) 


where W, = the fundamental frequency of the rotor shaft and rotor mass 
«11 = the fundamental frequency of the rotor shaft by itself 
«22 = the natural frequency of the rotor mass mounted on the rotor shaft in 
the absence of other masses 
a2? = the influence coefficient due to an unit load applied at the location of the 
rotor disk 
m = the mass of the rotor disk. 


Solving equation (7) for the first lateral natural frequency yields 


2 
Des Wi} (8) 
PG ar9mw}] 


For a simply supported beam of uniform mass, M, the fundamental frequency of the 


rotor shaft is 


209 


or -‘{ 2) (9) 


The influence coefficient can be determined from the static deflection equation for a 


simply supported beam evaluated at z =a. 


Pba 
499 = y(z =a) = pe -(L? a? -? | (10) 


The first lateral natural frequency was determined while varying the location of the 


rotor disk utilizing a computer program for MATLAB. 


D. DISCUSSION OF RESULTS 

The theoretical bounded range of Wy versus disk position are plotted along with the 
experimental values, see Figure 11. As seen from the plot, the experimentally measured 
Wn fell between the bounded region for values of rotor disk position between five and nine 
inches. As the rotor disk position was moved towards either bearing support, the 
experimental Wy began to fall below the bounded range. The percent deviation from 
theoretical values increased as the rotor disk was moved closer to the bearing supports. 

The reason for deviation from theoretical values can be attributed to the flexibility 
(stiffness) of the bearing supports and their effect of the rotor assembly. When the rotor 
disk is placed near the center of the shaft, the stiffness and damping of the bearing supports 
apparently have very little, if any, effect on the wy of the rotor assembly and may be 
modeled as pinned connections. As the rotor disk is placed closer to a bearing support, the 
effects of the bearing supports’ flexibilities play a larger role in the vibration characteristics 
of the rotor assembly. Therefore, the bearings can no longer be modeled as pinned 
connections but as springs (added flexibility), which would now account for the energies 


of the bearing supports. 
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Additionally. with these extra stiffness terms, the first mode shape of vibration will 
differ from the typical first bending mode associated with the pinned-pinned model. The 
mode shape is most likely a combination of bending, bouncing and rocking modes, with 
the bouncing and rocking mode having more of an effect as the disk is moved closer to the 
support. Therefore, using the static deflection as a mode shape is not representative of the 
actual mode shape when the disk is placed near the bearing supports, thus, increasing the 
theoretical error. Though, the shape function selected need not have the exact shape of the 
mode shape to obtain a good approximation, the accuracy of the theoretical values increases 
as the shape function used most resembles the actual mode shape [Ref. 4]. 

The theoretical bounded range was symmetric about the rotor shaft center, the 
experimental values were asymmetric. It can be concluded that the vibrational properties of 
the two bearings supports are different. Different stiffness values for the two bearing 


supports add to the complexity of the mode shape. 
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First Lateral Natural Frequency as a Function of Disk Position 
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Figure 11. Theoretical and Experimental First Lateral Natural 


Frequency as a Function of Rotor Disk Position. The Rotor Assembly 
is in a Simply Supported, Single Disk Configuration, Rotor Length of 
14.5 Inches, Disk with a Two Gram Mass Imbalance. 
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V. SYNCIRONOUS WHIRL DUE TO A MASS IMBALANCE 


It was desired to develop a Mechanical Engineering vibration laboratory for the 
phenomena of synchronous whirl due to a mass imbalance. The objective was to 
experimentally measure and observe the phenomena of synchronous whirl due to a mass 


imbalance during subcritical, resonant and supercritical operation of the rotor assembly. 


A. THEORETICAL BACKGROUND: SYNCHRONOUS WHIRL 

The behavior of rotor-bearing systems was first analyzed by Rankine in 1869. He 
incorrectly predicted that the operating speeds could never exceed the first cntical speed 
[Ref. 5]. The phenomena of synchronous whirl excited by a mass imbalance was later 
addressed by Jeffcott [Ref. 6]. Jeffcott explained that while an unbalanced shaft rotates 
about its longitudinal axis, the shaft will bend and, thus, whirl about the bearing axis (AB), 
see Figure 12. As the name implies, the whirl speed is equal to the rotational speed of the 
rotor. 

The amount of bending or whirl radius is dependent upon the center of mass offset 
(q), rotating speed (w) and its nearness to Wp. The steady-state solution of the differential 
equations of motion for this system yields equal shaft displacements in the x and y 


directions, resulting in circular whirl orbits of radius (R) (Ref. 7]. 


(12) 





7A 


Synchronous Whirl Condition. 








Ideal Circular 
Rotor Orbit 


a rs 
Figure 12. Synchronous Whirl and Synchronous Whirl Orbit. 


Another characteristic associated with synchronous whirl is a phase shift between the 
mass imbalance (a fixed rotor position) and the maximum shaft displacement. As the shaft 
whirls, the maximum displacement will begin to lag the mass imbalance. The amount of 
phase lag (A) also depends on the rotating shaft speed and its nearness to the first lateral 


natural frequency. [Ref. 7] 
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( 
: Ca) (13) 


Typical plots of whirl radius or shaft displacement and phase lag versus rotating shatt 


speed are shown in Figure [3 fora system with slight damping. 
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Figure 13. Theoretical Plots of the Whirl Radius (R) and the Maximum 
Displacement-Mass Imbalance Phase Lag (A) Versus Rotor Speed for a 
System with Slight Damping. 
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B. EXPERIMENTAL CONFIGURATION 
1. Rotor Assembly 
The rotor assembly was in a simply supported, single disk configuration. The 
bearing supports were positioned 14.5 inches apart (L=14.5) with the outboard bearing 
support flushed with the rotor shaft end. The rotor disk was centered between the bearing 
supports, 7.25 inches from either bearing support, with the zero degrees position hole 
aligned with the single notched disk. A mass imbalance of 1.6 grams was placed at the 
zero degrees position. With the mass imbalance aligned with the single notched disk, rotor 
position corresponds to the mass imbalance position. For initial balancing of the rotor 
assembly, a 0.25 gram mass was positioned at the 315 degrees position. 
2. Proximeter Connections 
An X-Y probe mount was placed 1.5 inches from the inboard bearing support. 
This probe mount position ensures that the vibrating rotor shaft does not contact the 
proximeter during operation. The connections for the proximeters are listed in Appendix 
Cc: 
3. Stroboscope 
A stroboscope is required for precise measurement of rotor speed throughout the 
experiment. Reflective tape was placed on the face of the rotor disk mass. 
4. Virtual Instrument 
With the NI AT-MIO-16f-5 Mutifunction I/O board configured for differential 
input as discussed in Appendix D, the virtual instrument, ROTOR.VI, was utilized with the 


default input parameters. See Appendix E for a description of ROTOR.VI. 
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C. DISCUSSION OF THE RESULTS 

The VI was ran for speeds ranging from 380 RPM up to 7000 RPM. Increments of 
speed varied from 300 RPM to 20 RPM, utilizing the smaller increments for speeds near 
the higher amplitude signals. The rotor is driven by the drive motor in a clockwise 
direction, considered the positive direction. 

For each run, the rotor speed and the phase shift between the maximum y- 
displacement and mass imbalance were recorded, as well as, sending the raw data to the 
computer hard drive. Appendix F contains a list of rotor speeds and corresponding data 
files for which the experiment was conducted. 

1. Maximum Rotor Displacement Versus Rotor Speed 

The maximum y-displacements versus rotor speed were plotted in Figure 14. 
Compared to the theoretical expectations of Figure 13, the results resemble a typical plot for 
a system with slight damping. The first lateral natural frequency determined from peak 
displacement for this experiment was 2740 RPM. 

Observation of the rotor orbit on the front panel of the VI showed that the orbit 
was not circular but elliptic during operation near Wy. Therefore, the x-displacements and 
the y-displacements were not equal as originally expected. Utilizing the raw data, 
maximum x-displacement versus rotor speed was also plotted. This x-displacement 
indicated a critical peak at 2900 RPM. A plot of the major radius of the elliptical orbit (R) 
versus the rotor speed showed these two critical peaks or a double critical condition. See 


Figure 15. 
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Figure 14. Experimental Results of the Maximum Y-Displacement 
Versus Rotor Speed. The Rotor Assembly is in a Simply Supported, 
Single Disk Configuration, Rotor Length of 14.5 Inches, Disk Located 
7.25 Inches from the Inboard Bearing Support (Centered) with a 1.6 
Gram Mass Imbalance. 
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Figure 15. Experimental Results of the Whirl Orbit Major Radius (R), 
Maximum Y-Displacement and Maximum X-Displacement Versus Rotor 
Speed. The Rotor Assembly is in a Simply Supported, Single Disk 
Configuration, Rotor Length of 14.5 Inches, Disk Located 7.25 Inches 
from the Inboard Bearing Support (Centered) with a 1.6 Gram Mass 
Imbalance. 
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2. Rotor Orbit Characteristics 
Utilizing the raw data sent to the computer hardrive, the x-direction and the y- 
direction voltages where converted to zero mean displacements by subtracting away the DC 
voltage component and then converting to mils. This allowed for detailed observation of 
the rotor orbit caused by synchronous whirl during subcritical, critical and supercritical 
speeds. As stated earlier, the first critical speed was determined to be 2740 RPM. 
a. Subcritical Speeds 

For rotor speeds between 300 and 1800 RPM, the rotor orbit was nearly 
circular, centered about the bearing axis. The radius of orbit increased in magnitude as 
rotor speed was increased. The direction in which the rotor traveled around these circular 
orbits was clockwise, thus undergoing forward whirl. 

As speed was increased from 1800 RPM, the rotor orbit changed from 
circular to elliptic with the major axis of the ellipse lying 60 degrees from the x-axis, see 
Figure 16. With increasing speeds, the orbit became more elliptical, that is, the major axis 
increased in magnitude and the minor axis decreased in magnitude. The angle of the major 
axis continued to be approximately 60 degrees from the x-axis. The rotor continued to 
experience forward whirl and traveled around these elliptical orbits in a clockwise direction. 

b. Critical Speed 

At the first critical speed (2740 RPM), determined by peak y-displacement, 
the orbit was elliptic, with the major axis at its maximum magnitude and at an angle of 60 
degrees from the x-axis. See Figure 17. The rotor continued to experience forward whirl. 

c. Supercritical Speeds 

The first rotor speed measured after the first critical speed was 2760 RPM. At 
this speed, the rotor orbit nearly resembled an straight line with an angle slightly less than 
that of the critical speed orbit. With continued small increments of speed, the orbits again 


became elliptic with the major axis slightly decreasing in magnitude and the minor axis 
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increasing in magnitude. The rotor reversed its direction of travel and now traveled around 
the elliptic orbits in a counter-clockwise direction. Therefore experiencing backward 
whirl, which means the precession is in the opposite sense than the shaft spin. Further 
increases of speeds up to the second cnitical peak, 2900 RPM, showed the elliptic orbits 
once again increasing its major axis and decreasing its minor axis until the orbit, once 
again, approximately resembled a straight line path. It is here, at the second critical peak, 
that the rotor again reverses its direction of orbital travel and reverts to traveling around the 
orbits in a clockwise direction. Once again undergoing forward whirl. This phenomena of 
two whirl direction reversals occurred within approximately 160 RPM following passage 
through the first critical speed. The angle of the major axis continued to rotate towards the 
x-axis. See Figure 18. 

As rotor speed was increased further, the orbits began to become less elliptic 
and continued to rotate their major axis towards the x-axis. This continued until the orbits 
became nearly circular once again at approximately 4000 RPM. See Figure 19. 

The orbits remained circular with a converging radius for the remaining runs 
with the exception of speeds near twice Wp. As seen from Figure 20, at 5560 RPM the 
rotor is experiencing synchronous and subsynchronous whirl. A power spectral density 
analysis of this signal indicates the synchronous frequency of 5560 RPM and the 
subsynchronous frequency of the first lateral natural frequency, Wp (2740 RPM). This 
phenomena follows the characteristics of an oil whip condition, also seen in the cascade 
plots described earlier in Chapter III. 

3. Maximum Y-Displacement - Mass Imbalance Phase Shift 
The phase shift between the y-displacement and the mass imbalance position was 
plotted versus rotor speed, see Figure 21. The experimental results resemble a typical plot 
for a system with slight damping. For rotor speeds <<wp (2740 RPM), the displacement 


and mass imbalance were in phase. As the rotor speed was increased towards wn, the 


os. 


displacement curve began to lag the mass imbalance. The phase shift increased slightly at 
first (30 degrees in 800 RPM's) but then increased sharply as the rotor speed approached 
Wn. At Wp, the phase shift was measured to be approximately 90 degrees lagging, as 
expected. As speed was increased further, the phase shift continue to increase sharply to 


160 degrees lagging and then slowly approached 180 degrees lagging. 


D. BACKWARD WHIRL 

The phenomena of backward whirl is explained by Vance using a simple two degree of 
freedom model with asymmetric, uncoupled bearing support stiffnesses. His model shows 
that the two critical speeds result from the uncoupled. asymmetric stiffnesses of the 
bearings. Backward whirl is experienced when the rotor speed is between the two critical 
peaks, as shown by the time rate of change of the whirl angle being negative or of opposite 
sense of the direction of shaft rotation. [Ref. 8] 

A more complex analytical model along with supporting experimental results was 
published by Subbiah, Bhat and Sankar. They used a six degree of freedom model with 
coupled, asymmetric bearings. Their experimental set-up contained a rotor with two 
identical hydrodynamic bearings at the ends. Again, they showed that the phenomena of 


backward occurred between the two critical speeds. [Ref. 9] 
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Rotor Orbits for 1809 - 2465 - 2556 - 2616 - 2675 - 2715 RPM 





1.5 Stolle econo! stiete\oie sles sis sislielieee Vee els sieioicie rei ieia suse piel a eee, «F515! sels ersie evel acelevevereierere 7 sie sac erehel sone eieuchele le elev nels evcilere ee 018 ses) 6 aie @ eevenese 
Zh is Ny 
a4 

UR unre ee ener are re co 7 ener ee es, See 2675 hh) CR to ee nee nee 

oO 0.5 abe ievalte felarenel olevenciecececeroner siete ovens eevone <r evckeerone stor snelssaNerierei oe cuter S tf ue VLG rete Te RST OMRON aes -dekg a Oia a vete Lose cerl’ 6f.ei oe: oie iene) sun tareterereuees! © 
= 
~— 
= 
2 

= (9) eaoeeceprenecre cer tocrr eer coreerrcetr orn) (ic eerie ee yy / Maret eerie eens meron 
O 
S 
ae 
L 
1s 

>»~-0.5 ii Wiese bolic oisia pW ain le atale cole este nig sve sialnievauee cwittesiie BUR cM cle ste © anid SECM TOE « ¥. 6 GAs 66 bis 34:4 6S: @ 9. #9 WN Be Gale lals S O'S Welw eig aielw © eaene 

BL Rear inn cur Secret 5° 6 | | Aer tay” A Se an eee ee ee ee ee nn ene ee eee ee 

1.5 eteoeoeveeveeeeee eee ee ees pe ese sae eo 8 8, Sie ie e808 Sah inde Sa iey Same eRe RES ce eoeeeeee eee eevee ee Busse 8 (exe vere re sicee exeseret eoeeeteeerereeet eee eee pone 

1.5 -4 -0.5 0 0:5 1 ine 


x-displacement (mils) 


Figure 16. Rotor Orbits for Speeds Subcritical Speeds of 1809-2465 - 
2556 - 2616 - 2675 - 2715 RPM. The Rotor is Experiencing Forward 
Whirl. The Rotor Assembly is in a Simply Supported, Single Disk 
Configuration, Rotor Length of 14.5 Inches, Disk Located 7.25 Inches 
from the Inboard Bearing Support (Centered) with a 1.6 Gram Mass 
Imbalance. 
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Figure 17. Rotor Orbit for the First Critical Speed 2740 RPM. The 
Rotor is Experiencing Forward Whirl. The Rotor Assembly is in a 
Simply Supported, Single Disk Configuration, Rotor Length of 14.5 
Inches, Disk Located 7.25 Inches from the Inboard Bearing Support 
(Centered) with a 1.6 Gram Mass Imbalance. 
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Rotor Orbits for 2760 - 2784 - 2809 - 2840 - 2898 RPM 
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Figure 18. Rotor Orbits for Supercritical Speeds of 2760 - 2784 - 2809 
- 2840 - 2898 RPM. The Rotor is Experiencing Backward Whirl. The 
Rotor Assembly is in a Simply Supported, Single Disk Configuration, 
Rotor Length of 14.5 Inches, Disk Located 7.25 Inches from_ the 
Inboard Bearing Support (Centered) with a 1.6 Gram Mass Imbalance. 
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Rotor Orbits for 2952 - 3022 - 3189 - 2501 - 4104 - 4603 RPM 
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Figure 19. Rotor Orbits for Supercritical Speeds of 2952 - 3022 - 3189 
- 3501 - 4104 - 4603 RPM. The Rotor is Experiencing Forward Whirl. 
The Rotor Assembly is in a Simply Supported, Single Disk 
Configuration, Rotor Length of 14.5 Inches, Disk Located 7.25 Inches 
from the Inboard Bearing Support (Centered) with a 1.6 Gram Mass 
Imbalance. 
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Rotor Orbit - Rotor Speed 5560 RPM 
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Figure 20. Rotor Orbit for Supercritical Speed of 5560 RPM. The 
Rotor is Experiencing Synchronous and Subsynchronous' Forward 
Whirl. The Rotor Assembly is in a Simply Supported, Single Disk 
Configuration, Rotor Length of 14.5 Inches, Disk Located 7.25 Inches 
from the Inboard Bearing Support (Centered) with a 1.6 Gram Mass 
Imbalance. 
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Phase Shift vs. Rotor Speed 


50 





0 1000 2000 3000 4000 5000 6000 7000 8000 
Rotor Speed (RPM) 


Figure 21. Maximum Y-Displacement - Mass Imbalance Phase Shift 
Versus Rotor Speed. The Rotor Assembly is in a Simply Supported, 
Single Disk Configuration, Rotor Length of 14.5 Inches, Disk Located 
7.25 Inches from the Inboard Bearing Support (Centered) with a 1.6 
Gram Mass Imbalance. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

An experimental facility was developed to investigate the rotordynamic phenomena of 
rotating machinery. The facility consists of a Bentley Nevada Corporation Rotor Kit with a 
computerized data acquisition and control process. A LabVIEW Virtual Instrument was 
designed to control the data acquisition of the rotor position and rotor displacements. 
Experiments conducted included construction of Cascade plots, analytic prediction of the 
first lateral natural frequency while varying rotor disk placement and demonstration of a 
synchronous whirl condition during which backward whirl was experienced. 

As seen from the results of the various experiments conducted, it can be concluded that 
the experimental facility can effectively demonstrate fundamental rotordynamic phenomena 
during subcritical, resonant and supercritical operation and serves as test for such. The 
experimental facility allows for the determination of the rotordynamic characteristics 
associated with various rotor configurations, as well as, the observation and analysis of 
unique rotordynamic phenomena as configuration parameters are changed. 

The Nevada Bentley Rotor Kit is a versatile assembly and can easily be configured to 
model many phenomena that occur in rotating machinery. The flexibility of Lab VIEW's 
virtual instruments allow for precise data acquisition control that can easily be applied to 
each rotor assembly configuration. The ability to collect, process and store data, enabling 
detailed analysis, has been demonstrated. 

In summary, in consideration of the experiments conducted and the analysis of the 
data, it can be concluded that the experimental facility can be utilized to enhance 
rotordynamic research, as well as, serve as a valuable educational tool. Utilizing an model 


of a rotordynamic system that effectively duplicates the rotordynamic phenomena of its 
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larger actual system, allows for easier, more cost effective rotordynamic analysis. This can 
ultimately result in the prediction, elimination and/or avoidance of undesirable, machine 
damaging rotordynamic phenomena 

Other findings that can be drawn from the supporting experimental data are directed at 
the rotordynamic characteristics of the simply supported, single disk rotor configuration 
utilized for each experiment. In particular, the facility demonstrated: 

I. The construction of Cascade plots allow for the determination of general 
fundamental rotor characteristics and the rotordynamic phenomena experienced during 
subcritical, resonant and supercritical operation of the rotor assembly. Synchronous 
response and other subharmonic or superharmonic responses caused by nonlinearities can 
be detected, observed and analyzed from this plots. 

2. The rotor assembly can be modeled as a pinned-pinned connection when the 
rotor disk is located near the center of the shaft. The experimental first lateral natural 
frequency can be analytically bounded using Rayleigh's Energy Method for the upper 
bound and Dunkerley's method for the lower bound when configured in this way. 

3. A pinned-pinned analytical model of the rotor assembly is not valid when the 
rotor disk is located near the bearing supports. The stiffness effects of the bearing supports 
must be accounted for in modeling. Additionally, the stiffness properties differ for each 
bearing support. 

4. The unique rotordynamic phenomena of backward whirl was demonstrated 
during the synchronous whirl experiment. Cause of this phenomena was not investigated 
but was assumed to be caused by the cross-coupled or asymmetric properties of the 
bearings. This could be attributed to oil soaking of the outboard bearing in standard 
automobile grade lubricating oil prior to running the experiment. Additionally, the effects 
of a mass imbalance can easily be seen in this experiment, thereby, reinforcing the 


importance of machine balancing. 
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5. The rotor assembly experienced oil whirl and oil whip conditions. This is 
evident in both the cascade plots and the synchronous whirl experiment. It should be noted 
that the experimental runs for the cascade plots were conducted prior to soaking the 


outboard bearing in the oil and, still, the rotor assembly experienced oi] whirl and oil whip. 


B. RECOMMENDATIONS FOR FUTURE WORK 

This research has identified several areas that warrant further studies. The most 
obvious is to continue to determine the rotordynamic characteristics of the Nevada Bentley 
Rotor Kit using different rotor assembly configuration. Examples of various rotor 
assembly configurations included in the kit are overhangs, oil whirl and free spinner 
options, double disk configurations and rub conditions. 

Additionally, an important compliment to the experimental facility would be an precise 
analytical model of the rotor assembly. The model should not only accurately predict the 
natural frequencies and mode shapes of the rotor assembly for various rotor configurations 
but also predict the occurrence and rotor response of such rotordynamic phenomena as oil 
whirl, oil whip, and backward whirl. To build an effective analytical model. the bearing 
characteristics, as well as, other assembly components should be experimentally 
determined. 

The fluid-structure interaction problem can be investigated experimentally in the 
facility. Accurately determining the properties of the Oilite bearings and the effect of oil 
soaking the bearings or eliminating oil whip and oil whirl conditions are just a few of the 
applicable topics. 

Overall, the experimental facility can be effectively utilized for continued rotordynamic 


research, as well as, educational purposes. 
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APPENDIX A. LIST OF DIMENSION AND PROPERTIES 


Rotor Shaft: 


Material: 4140 Stainless Steel 

Modulus of Elasticity: 200 GPa 29.2(10) psi 

Length: 45.7 cm 18.0 in 

Diameter acorn 0.375 in 

Mass 0.253 kg 0.558 Ibm 

Density 7766 ke/m3 0.281 Ibm/in3 
Rotor Disk: 

Material: 316 Stainless Steel 

Width: 2.54 cm 1.0 in 

Inner Diameter 0.9525 cm 0.375 in 

Outer Diameter ) thes) ie 3.0 in 

Mass 0.816 kg 1.8 Ibm 
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Proximeter Calibration Curve 


APPENDIX. B 
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Figure 22. 
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APPENDIX C. PROXIMETER CONFIGURATION 


For all experiments conducted, the electrical connections remained the same. A listing 


of the electrical connections for each proximeter is found in Table IT. 


A. EDDY CURRENT PROXIMETERS 

An X-Y probe mount was utilized for all experiments. Location of the probe mount 
must ensure that the vibrating rotor shaft does not contact the proximeters during operation. 
Prior to operation, each proximeter must be positioned with a separation distance from the 
rotor shaft of approximately 40 mils. This can be accomplished by disconnecting the motor 
power supply lead, turning the power supply on and reading between -7 and -8 volts DC at 
the proximeter. 

Severe inference results when using an X-Y probe stand and placing the two probes in 
the same plane. No interference was found when using arn X-Y probe stand with the two 
proximeters in different planes. 

It should be noted that the output signal from the twelve notch disk proximeter is 
connected directly to speed controller and not connected to a proximeter transducer like the 
other proximeter outputs. This is due to the fact that the twelve notched disk provides 
feedback for the motor speed controller. This signal can be acquired from the proximeter 
(PROX) and common (COM) terminals on the back of the speed controller. 


B. TTL SIGNAL CONVERTER CONNECTIONS 

The TTL Signal Converter converts the single notched and 12 notched disks 
proximeter outputs to a TTL signal. The TTL signal converter's power supply lead (red) is 
connected to the +5V terminal on the SC-2070. The TTL converter's ground lead (green) 
is connected to the DGND terminal on the SC-2070. 
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TABLE I. ELECTRICAL CONFIGURATION 


Proximeter Proximeter TTL SC-2070 Channel . 
Probe Transducer | Converter differential (single) 


l (1 - DGND) 
Single notched l white in (8 - white out) 


y (2 - common) 
X - direction (9 - output) 


Se a. “(0 eutgu 
- direction = (10 - output) 
frveenowhed | nla tr back 
Twelve notched n/a (11 - black out) 
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APPENDIX D. NI AT-MIO-16F-5 CONFIGURATION 


Following installation ot the NI AT-MIO-16F-5 mutifunction I/O board and 
installation of the NI LabVIEW software package, software configuration of the NI AT- 
MIO-16F-5 was required prior to operation. LabVIEW needs the software information to 
recognize the NI AT-MIO-16F-5 and to set default data acquisition parameters. 

Configuration is performed through the use of NI-Data Acquisition (NI-DAQ) 
Configuration Utility version 4.4.2. This utility is accessed through the WDAQConf icon 
in the LabVIEW Program Group of the Windows Program Manager window. Detailed 
instructions are found in [Ref. 10}. 

The NI AT-MIO-16F-5 is configured as shown in Table HI with default data 


acquisition parameters listed in Table IV. 


TABLE III. NI AT-MIO-16F-5 CONFIGURATION 


TABLE IV. NI AT-MIO-16F-5 DEFAULT DATA ACQUISITION 
PARAMETERS 


DAQ Mode Analog Input 
Polaritv/Range 0 to +10 volts 
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APPENDIX E. ROTOR.VI VIRTUAL INSTRUMENT 


A. INPUTS 
With the NI AT-MIO-16f-5 Mutifunction I/O board configured for differential input as 

discussed in Appendix D, the following is a description and listing of the front panel inputs 
required for execution of the VI. Default values are listed in parenthesis. See Figure 22 for 
a copy of the front panel. 

1. CHANNELS: The analog input channels. (1, 2, 3, 4) 

2. NUMBER OF SCANS (N): The number scans to acquire for each channel. (2048) 

3. SCAN RATE (scans/sec): The sampling frequency, the number of scans to acquire 

per second. (4096) 

4. INPUT LIMITS (volts): Allows data to be acquired that is outside the polarity 

settings of the AT-MIO-16F-5 board. The board is configured for 0-10 volts. 

(high limit: 15, low limit: Q). 

5. TRIG CHANNEL: The channel to be used as the trigger. (1) 

6. PRETRIGGER SCANS: The number of pretrigger scans to acquire. (0) 

7. TRIG SLOPE: Select triggering based on a rising or falling slope. (1: nsing) 

8. TRIG LEVEL (volts): The triggering condition, the voltage level that must be 


crossed for triggering to occur. (0.8) 


B. EXECUTION 
To execute this VI, click on the run button on the run mode palette with the mouse. 
When executed, this VI starts acquiring data from the specified analog input channels into a 


circular buffer. The size of the buffer is determined by the VI as either 10K scans, 0.4 


Dy | 


seconds or 2.5 times the number of scans to acquire, whichever is greatest. As data is 
acquired, it begins to search through the data for the trigger condition. Once the tngger 
condition is met, the retneval of data begins at the inputted scan rate. The data acquisition 
sub-virtual instrument (SubVI), called AI Read, acquires the data in a two dimensional 
matrix, with the channels in columns and the scan number in rows. Data acquisition 
continues until the specified amount of data is obtained. The data processing portion of the 
VI then begins. If the trigger condition is not met within a time limit of two seconds, the 
VI stops the execution of the VI and displays a TIMED OUT error measure on the front 
panel. 

The VI's block diagram and hierarchy display of the SubVI's utilized are found in 
Figures 23 and 24 respectively. 

To obtain a print out of the front panel following execution of the VI, ensure the print 
mode button on the run mode palette is active prior to running the VI, click on print icon 
with mouse to activate. Additionally, ensure the printer is set up for Black/Gray Scale 


Only. 


C. DATA PROCESSING, DISPLAYS AND OUTPUT PATIIS 

Essentially, the VI allows the data to flow two paths. One being a data processing and 
display path and the other an output file path . 

a. Data Processing and Display Path 

In the data processing and display path, the VI converts a copy of the raw data 

from voltage to mils for the acquired single notch, x-displacement and y-displacement data 
(channels 1, 2, and 3). Additionally, the VI determines and subtracts the mean for the x- 
displacement and y-displacement data (channels 2 and 3). This processed data is displayed 


on the front panel as a Rotor Orbit Plot and a Time Signal Trace. 


DIP 


The Rotor Orbit Plot is a graph of the y-displacement versus the x-displacement 
channel. Since the X-Y probe stand provides two measurements spaced 90 degrees 
circumferentially, rotor orbits or Lissajous patterns can be obtained. It is the path travel by 
the rotor during operation. Characteristics of the rotor orbit over a range of speeds can be 
observed from this plot. The limits for each axis are locked for this plot. 

The Time Signal Trace is a voltage versus time graph for the single notched 
channel and the y-displacement channel. Maximum y-displacement data can be obtained by 
moving the y-Displacement cursor with the mouse until aligned with peaks of the y- 
displacement curve. Movement of the cursors can be achieved by placing the mouse 
directly on the cursor and dragging the mouse or by clicking on the cursor direction arrows 
below the plot. The values of the cursor crasher position are displayed below the plot on 
the cursor displays. 

Additionally, phase shift between the y-displacement and single notch (mass 
imbalance position) traces can be obtained for each run of this VI. Selection of the y- 
displacement instead of the x-displacement allows for easier phase shift measurements; 
since the y-displacement leads the x-displacement by 90 degrees, the y-displacement begins 
a new cycle at the position of the single notched disk when in phase. 

As seen from the front panel display, Figure 25, the single notched trace is 
superimposed with the twelve notched trace. This is a result of signal interference that 
occurs in the TTL converter. Nevertheless, this allows for increments of 30 degrees on the 
single notched trace. 

The limits for each axis for the Signal Time Trace Plot are not locked and may be 
changed when necessary. This is accomplished by highlighting the desired limit to be 
changed with the mouse and inputting the new limit via the key board. The plot will 
instantaneously be replotted with the new limits when the mouse is click away from the 


inputted limit. 


yes: 


B. OUTPUT PATH 

For the output file path, the VI formats a copy of the acquired raw data to a spread 
sheet format with the first dimension (rows) being the scan number and the second 
dimension (columns) being the channel. This spreadsheet formatted data, channels 1 thru 


4, is then sent to an user specified ASCII file. 
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eee tng channel (1) R | 
c —— otor Orbit 
2.3, Baas Y vs. X Dispiacement (mils) 
number of scans scans {0} 7 
to acquire a0 | 
wl2024_ | a 
trig slope 
scan rate ae 7) 4] 
(scans/sec] 
=4(4096.00 _| 


input fimits (volts) 
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Figure 23. ROTOR.VI Front Panel. 
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Figure 24. ROTOR.VI Block Diagram. 
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Figure 25. ROTOR.VI Hierarchy. 


APPENDIX FF.) SYNCHRONOUS WITIRL TEST MATRIX 


Byte: 


File Name File Name 
Rotor Speed Rotor Speed 

C:\LabVIEW\ (RPM) C:\LabVIEW\ (RPM) 
MEDATA\ MEDAT A\ 
LABI1 380 LAB30 2809 
LAB2 53 LAB31 2840 
LAB3 706 LAB32 2898 
LAB4 887 LAB33 2952 
LABS O25 LAB34 3022 
LAB6 i ae | LAB35 3095 
LAB7 jk UY LAB36 3189 
LAB8& 1409 LAB37 3266 
LAB9 1519 LAB38 336m 
LAB10 1642 LAB39 3440 
LAB1I1 1809 LAB40 3501 
LAB12 1886 LAB41 36238 
LAB13 1998 LAB42 3800 
LABI14 2114 LAB43 3933 
LABI5 2207 LAB44 4104 
LAB16 2270 LAB45 4307 
LABI17 2326 LAB46 4425 
LAB18 2407 LAB47 4603 
LABI19 2465 LAB48 4802 
LAB20 E12 LAB49 4964 
LAB21 2556 LABSO 5 ie 
LAB22 2361 6 LABS 1 532 
LAB23 7105) LABS52 5360 
LAB24 Z2OTS LABS53 5794 
LAB25 2694 LABS4 6009 
LAB26 ge 15 LABS55 63:20 
LAB27 2740 LABS56 6574 
LAB28 2760 LABS7 6809 
LAB29 2784 LABS8& 7000 


(Ref. 1] 


{Ref. 2] 


[Ref. 3] 


[Ref. 4] 


[Ref. 5] 


[Ref. 6] 


(Ref. 7] 


[Ref. 8] 


{Ref. 9] 


(Ref. 10] 
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